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, Both linear and nonlinear Langevin equations are derived directly from the Liouville equation 

' for an exactly solvable model consisting of a Brownian particle of mass AI interacting with ideal 

gas molecules of mass m via a quadratic repulsive potential. Explicit microscopic expressions for 
(j^ ', all kinetic coefficients appearing in these equations are presented. It is shown that the range of 

' applicability of the Langevin equation, as well as statistical properties of random force, may depend 

not only on the mass ratio m/M but also by the parameter Nm/M, involving the average number 
A'" of molecules in the interaction zone around the particle. For the case of a short-ranged potential, 
when N <^ 1, analysis of the Langevin equations yields previously obtained results for a hard- 
wall potential in which only binary collisions are considered. For the finite-ranged potential, when 
• multiple collisions are important (A*' S> 1), the model describes nontrivial dynamics on time scales 

that are on the order of the collision time, a regime that is usually beyond the scope of more 
phenomenological models. 
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. I. INTRODUCTION 

The explicit derivation from first principles of the Langevin equation, describing the evolution of a small number 
of variables in a complex system, is often necessary since in many cases the statistical properties of the random 
force, the range of applicability and even the form of the equation are far from evident [1] . The need for microscopic 
[ considerations is especially compelling in the case of a nonlinear Langevin equation where the usual phenomenological 
(~| approach of adding stochastic terms to the deterministic nonlinear equation describing the relaxation of targeted 
O variables may be inadequate [2]. Although statistical properties of the random force in the nonlinear Langevin 
^ O ^ . equation may be deduced phenomenologically in some cases [3] , it is generally necessary to start from the microscopic 
behavior of the system in order to construct the appropriate form of the equation. 
, The conventional systematic method for deriving the Langevin equation (LE) for a Brownian particle exploits 
^ > Mori's projection-operator techniques [1], which allow the transformation of the microscopic Liouville equation to 
' a non-Markovian predecessor of the LE, generally known as the generalized Langevin equation. The LE can be 
J"' . then obtained by a subsequent perturbation expansion of the memory kernel appearing in the generalized Langevin 
' equation using the square root ratio of the mass of a bath particle to that of the Brownian particle, (m/M)-'^/^, as a 
, perturbation parameter A. While the first step in the derivation of the LE, involving rearrangement of the Liouville 
' equation with projection operator methods, is an exact algebraic procedure, the question of convergence of the A 
, expansion in the second step is subtle and can be strictly justified only under assumptions which are difficult to 
"Y^ • prove in general [4]. Mazur and Oppenheim developed an alternative projection-operator approach more suitable to 
' analyse the convergence of the A expansion [5] . The validity of the perturbation analysis in both the Mori and the 
Mazur-Oppenheim approaches has been examined for an exactly solvable model consisting of tagged particle motion 
in a harmonic lattice [6] . Unfortunately, the convergence properties of the A series are trivial for this model, since 
all terms in the A expansion higher than zeroth order vanish. As a result, the memory kernel for the linear damping 
force in the LE does not depend on the mass ratio, and all non-linear terms are identically zero [7]. 
Q One goal of this paper is to present and analyze a simplified model that serves as a useful and non-trivial testing 

I!* ground for examining some of the subtle aspects of the theory of Brownian motion. The model considered here 
, is a generalization of the well-known Rayleigh model of a Brownian particle of mass M constrained to move in 
■ one dimension and subjected to collisions with an equilibrium ideal (non-interactive) gas of particles of mass m. The 
5h ' Rayleigh model is perhaps the oldest model of non-equilibrium statistical physics [8] , and has attracted much attention 
over the years, with early work [9-11] focusing on the model as a test of the systematic derivation of macroscopic kinetic 
equations for the heavy particle from the master equation. More recent investigations have examined the stationary 
and transient solutions of the asymmetric Rayleigh model in which the thermodynamic parameters characterizing the 
gas to the left and right of the piston differ [12]. In all these studies, the interaction between the Brownian particle 
and the bath is assumed to be short-ranged with a negligible collision time Tc- Only binary collisions are considered 
in this model because the range of interaction is assumed to be short compared to the average distant between bath 
particles. The usual starting point for analysis of the Rayleigh model with binary collisions is a Markovian master 
equation for particle's velocity distribution function. The master equation is only an approximate form of the fully 
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microscopic Liouville equation, and valid only for time scales longer than Tc. To test many aspects of the theory 
of Brownian motion, it is essential to start from a fully microscopic description of the dynamics directly from the 
Liouville equation so that any non-Markovian character of bath correlations is properly incorporated and particle 
dynamics on time scales less than Tc can be described. To this end, one may generalize the interaction between the 
bath and Brownian particles from a hard- wall to a parabolic repulsive potential. The generalized model is analytically 
solvable, allowing explicit calculation of all terms appearing in the derivation of both the linear and nonlinear LE 
beyond the binary collision approximation. We demonstrate that in addition to the small parameter A = (m/M)^/^, 
the character of the dynamics of a tagged particle can be governed by an additional parameter, iV^/^A, involving the 
average number N of bath particles simultaneously interacting with the Brownian particle. For a large Brownian 
particle, it is shown that when N 1 and multiple collisions are important, the parameter of the formal A-expansion 
is actually N^/'^X. 

The paper is organized as follows: In section 2, the Mazur-Oppenheim approach is reviewed to provide groundwork 
for all subsequent analysis. In section 3, the structure of terms in the A expansion is examined and presented in 
a convenient form. In section 4, the general formalism is applied to the Rayleigh model with a quadratic repulsive 
potential describing bath-Brownian particle interactions, and the LE is derived for the heavy particle. The nonlinear 
LE is obtained in section 5 and various aspects of this equation are discussed. Finally, a few concluding remarks are 
made in section 6. 



II. BASIC EQUATIONS 

The Hamiltonian for a Brownian system composed of a tagged particle of mass M in a bath of point particles of 
mass m is 

H=^+H., (1) 

i 

where x = {xi} and pi are positions and momenta of bath particles, X and P are those of the Brownian (or tagged) 
particle, and Hq is the Hamiltonian of the bath in the field of the tagged particle fixed at X. One can expect that on 
average P ~ y/MhsT, where ks is Boltzmann's constant and T is the temperature, and that the scaled momentum 
P* = AP, where A = y/m/M, is of the same order as the average momentum of a bath particle. In terms of scaled 
momentum, one can write the Liouville operator as 

£ = £o + A£i, (3) 

where Fi = —VxiU and F = — VxU are the forces on the i-th bath particle and on the Brownian particle, respectively. 
The operator Cq dictates the dynamics of the bath in the field of the fixed Brownian particle. 

If the mass of the tagged particle is large (i.e. a Brownian particle), one might intuitively expect that inertial effects 
of the particle's motion are small and that the force on the particle, F{t) = e^*P, is close to the pressure force, i.e. 
to the force on the fixed tagged particle, 

Fo{t) = e^o'F. (6) 

In the Mazur-Oppenheim approach [5], the force F{t) is decomposed using the projection operator V which averages 
a dynamical variable A over the canonical distribution p = Z~^exp{—pHo), for bath variables at fixed position of the 
tagged particle, 

VA = {A) = J P^H dxi dpi, (7) 

i 

where Z is the canonical partition function and /S = l/ksT. Using the operator identity [13] 

^iA+B)t ^ ^At ^ /•* dre^(*-)Se(-^+^)^ (8) 
Jo 
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with A = C and B = —VC, one may formally decompose the force F{t) on the tagged particle into a "random" part 
and a remainder as 



Jo 



(9) 



where FT{t) = e^^'F and Q = 1 - P. The factor VCF^t) in the integral in Eq. (9) can be simplified taking into 
account the orthogonality of V and jCq (i.e. VjCq = 0), and the equality 

{VxFHt)) = -piFoF^t)), (10) 

which can be derived by integration by parts. As a result, one obtains the following exact equation of motion for the 
scaled momentum of the tagged particle, 

^ = XF^t) + 1^ dr e^(*-) (Vp. - {FF^r)), (11) 

where F\t) is a zero-centered random force obeying (F^(t)) = Ve^^*F = 0. 

The random force F^{t) = e^'~°~^^^'~^^^F can be further expanded in terms of the mass ratio parameter A using 
identity (8) to obtain 

F\t) = Fo{t) + A / dh e^°(*-*i)Q£iFo(ti) + (12) 
Jo 

A^ / dh / '(^^2e^°(*-'l)Q£le^°('^-'=)Q£li^o(^2) + ••• 
Jo Jo 

Similarly, the kernel K{t) = {FF^{t)) appearing in the exact equation of motion (11) of the tagged particle may be 
expanded in a power series in A 

K{t) = {FF^{t))=^Ki{t), (13) 
I 

Ko{t) = X''{FFo{t)), 

K^{t)=X^ f dtiCi{t,ti), 
Jo 

K2{t) = X^ [ dti I dt2C2{t,ti,t2)--- 

Jo Jo 

where the correlation functions C'l are defined to be 

Ci{to,h) = (i^ (e^°(*«-*^)Q£i) Fo{h)) , (14) 
C2{to,ti,t2) = {f (e^o(*''-*i)Q£i) (e^''(*i-*^)Q£i) Fo{t2)] 

Ci{to,hr--,ti)=/F(jJe''''^'--'-'-^QCi)Fo{ti)\. 



i=l 

The truncation of the A-expansion to zeroth order, K{t) w Ko{t), leads from Eq.(ll) directly to the generalized 
Langevin equation 

= AFt(i) - X' dr Mo{T)P4t - r), (15) 

where 

Mo{t) = Ixoit) = j^iFFoit)). (16) 

This approximation is sensible provided the correlation functions C; appearing at higher order in the A-expansion 
decay on a similar A-independent time scale Tc characteristic of motions of the fixed-particle system (i.e. governed 
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by Co). Mazur and Oppenhcim [5] succeeded in proving that this is the case assuming the factorization properties 

{A{t,)e^otB{t,)} {A{t,)){B{t,)). 

While being formally non-local in time, Eq. (15) can actually be written in a form that is local in time by expanding 
P* [t — t) around r = to obtain 

^=XF\t)-\'^o{t)P.{t\ (17) 
where 70 (t) = /g drMoij). The non-local correction terms to this approximation are of the form [14] 

\^ ( drMoir) I dr' P^t') X^, (18) 

^0 Jt-T 



and are therefore of higher order in A. Naturally, this analysis is pertinent only if the characteristic time Tc for decay 
of Mo{t) docs not depend on A. 

The local equation (17) is applicable on arbitrary time scales and for t ^ Tc assumes the form of the conventional 
LE with a time- independent damping coefficient 70 = dtMo{t). It is then evident from (17) that the autocorre- 
lation function of the momentum of a heavy particle decays on a time-scale Tp ^ A~^ that is much longer than the 
characteristic time of the bath Tc ~ A". One can therefore expect the local form of the LE with time- independent 
damping to be a good approximation for Eq. (17) except at short times determined hy t < Tc- 

It will be shown below that for homogeneous bath Ki{t) = 0, so the next approximation for the A-expansion (13) 
is of the form K{t) ~ Ko{t) + K2{t). The equation of motion for the momentum of the tagged particle in this case 
includes a non- linear damping term of third order in P*. This equation will be considered in section 5. 



III. STRUCTURE OF THE TERMS IN A-EXPANSION 



To examine the convergence properties of the A-expansion and other features of the projection operator derivation 
of the LE, the structure of the correlation fimctions Ci defined in (14) which appear in the expansion of the memory 
function must be analyzed. Although only the functions Cq and C2 are needed to obtain the nonlinear LE to lowest 
order in A, it is useful to know general properties of Cj. 

By inspection of the symmetry properties of the system, it is immediately apparent that the correlation fimctions 
C^n+i corresponding to odd powers of A contain an odd number of factors F and Vx, and therefore vanish for isotropic 
systems. In fact, in the absence of external field, the dependence on the particle coordinate appears only through the 
difference Xi — X, and it is useful to introduce new variables qi= Xi — X. Since the vectors F = — WxU = J^i ^qi^ 
and Vx = — J2i ^qi have negative parity and the Hamiltonian Hq is invariant with respect to transformation 
{qi —Qi, Pi —Pi}, the correlation functions C2n+i{t,ti, ...t2n+i) vanish. 

Correlation functions of even orders C2n do not vanish and have a rather complicated structure. For notational 
simplicity, we restrict the analysis to the case of one-dimensional diffusion, expecting no physical features in higher 
dimensions. For future development, it is convenient to define 

Go{t) = Fo{t), 

Giith) = S{t-h)F„ih), 

G2{t, h,t2) = Sit - h)Sih - t2)Fo{t2), 

G,(i, ti, • • • , ts) = S{t - ti) • • • S{t,-i - t,)Fo{ts), (19) 

where 

5(t,-fO = e^''(*^-*''^. (20) 

Note the property 

{Gi,Gi,---Gi,)=0 (21) 

which holds for arbitrary time arguments when ii + i2 + ■ ■ ■ + ii and I have different parities. For example, (Gj) and 
{GoGoGi) vanish for even i, while (GoGj) and (GoGoGoGi) are zero for odd i. 
Using the definitions above and according to Eqs. (14), the second order correlation function C2 can be written as 



4 




(22) 



Using cumulants, denoted by (( • • • )) and defined via the relations, 



{A1A2) 



((^1)), 

{A^){A2) + iA^A2)), 
{A,){A2){A^) + 

{A{){{A2A^)) + (^2)((^3^l)) + (^3)((^2^l)) + 

{{A^A^As)), 



the function C2 can be written as 



C2{tM,t2) = (^-^j {{GoG2{t,tl,t2)}} + 

-{{GoGo{t-h)Gi{t,t2))). 




(23) 



Note that the zeroth order kernel Kq is also a cumulant, Kq = {FFo{t)) = ((GoGo(t))). The relevance of cumulant 
representation for C; follows from the fact that one expects the cumulants to have similar scaling properties with 
respect to parameters of the system independent of their order. For example, it will be established in the next 
section that cumulants {{Gi^ Gi^ ■ ■ - Gi^ )) of any order are linear functions of the average number N of particles in the 
interaction zone around the particle for the Raylcigh model. Therefore, for this model, the first two non-vanishing 
terms in the A-expansion (13), namely Kq and K2, are both linear in N. 

Consider the correlation functions C21 with 1>1. It is tedious, though not difficult to establish that these correlation 
functions are of order I in cumulants {{Gi^Gi^ • • • Gj,)). Note that C2i has the contributions which are of order I + 2 
or less in Gj. First consider the terms of maximal order 1 + 2: One type of these terms are the correlation functions 
of the form 



where here and below time arguments have been omitted for brevity. The indices ii, - ■ ■ ,ii may take any values from 
the set {0, 1, • • • , 1} provided 



For example, for I ~ 1 the correlation function of maximal order in Gj which contributes to G2 is (GoGqGi), as can 
be seen from Eq. (23). For / = 2 the function G4 includes contributions of the form (G0G0G0G2) and (GoGoGiGi), 
and so on. 

All other contributions of order 1 + 2 in Gi can be written as products of correlation functions of lower orders which 
can be obtained dividing the sequences 



into groups in all possible ways without permuting functions. For example, for Z = 2 the three remaining terms are 



Clearly, the terms of maximal order in cumulants corresponds to the case when all indices ig in Eqs. (24) and (25) 
are equal to one. Then the cumulant expansion of all such terms contains the contribution 



{GoGoGi^Gi^ ■ ■ ■ Gil), 



(24) 



ii+i2-\ =1- 



(25) 



Go(ri)Go(r2)Gi, (r3)Gi, (T4) ■■■Gi, {n+2) 



(26) 



(Go(ri)Go(T2))(Go(r3)G2(r4)), 

(Go(ti)Go(t2))(Gi(t3)Gi(t4)), 
(Go(ti)Go(t2)Gi(t3))(Gi(t4)), 

(Go(Ti)Go(T2))(Gi(T3))(Gi(r4)) 
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(GoGo)(Gi)(Gi)---(Gi), 



(27) 



which contains I factors {Gi{ti)) and, therefore, is of order Z + 1 in cumulants. However, all terms of order Z + 2 in Gj 
described above enter in C2; via combinations 

{GoGoGi-^Gi^ ■ ■ ■ Gi^) — {GoGo){Gi-^Gi2 ■ ■ ■ Gi^}, 
(GoG'oGij)(Gi2 • • • Gi,) — (GaGo) {Gi-^) {Gi^ ■ ■ ■ Gi^), 
{GoGnGi-^Gi.^){Gi^ ■ ■ ■ Gi^) — {GqGq) {Gi-^^Gi^) {Gi^^ ■ ■ ■ Gi^), ■ ■ ■ 

which is a consequence of a presence of the operator Q (first from the left) in the definition of Gj (see Eq. (14)). 
As a result, the contributions (27) involving I + 1 cumulant factors cancel, and the maximal order in cumulants of 
surviving terms in the expression for G21 contain at most I cumulant factors. 

Having established that the contributions to G21 from terms of maximal order I + 2 in G, are of order / in cumulants, 
it is clear that the terms of order I + 1 and lower in Gj cannot contain more that I cumulant factors since such terms 
contain at least one isolated factor of Go whose average vanishes. Hence C21 contains at most I cumulant factors. 

In the next section, we examine the consequences of the cumulant expansions of the memory function for a specific 
system, namely, the Rayleigh model with a repulsive parabolic potential. It will be demonstrated that for the Rayleigh 
model, cumulants are linear functions of the average number A'' of particles in the interaction zone around the particle. 
This, in turn, implies that the terms in the A-expansion behave as 

{N + N'^ + --- + N%^^ (28) 

for ? > 1. Clearly, for a large particle and/or long-ranged potential leading to :» 1, these results suggest that 
K21 ~ A/^'A^', demonstrating that the actual parameter of A-expansion for this exactly solvable model is in fact ArV2A. 
On the other hand, for a short-ranged potential, when A ^ 1, all terms in the expansion arc linear in A^. and one 
sees that K21 ~ A^A^', and the effective small parameter of A-expansion is in fact the square-root of the mass ratio. 
To further illustrate this analysis, the explicit form of the cumulant expansion of G4 is presented in Appendix A. 



IV. THE LINEAR LANGEVIN EQUATION FOR A HEAVY PARTICLE OR IDEALIZED PISTON 

Consider the random motion of a piston of mass M and cross-sectional area S siibjcctcd to collisions with an ideal 
gas particles of mass to. The gas particles and the piston are constrained to move in one dimension perpendicular to 
the piston faces. The velocity distribution of incident particles fuiv) before collision with the piston is Maxwellian 
with inverse temperature /3, namely 

fM{v) = ( exp {-UrmA . (29) 




. 27r y V 2 

The piston-particle interaction is assumed to be described by a purely repulsive parabolic potential. For particles to 
the left of the piston the interaction potential between a gas particle and the piston is 

(30) 

where kf is a force constant, x is the coordinate of the gas particle, Xi = Xif — a the boundary of the piston-particle 
interaction zone, Xif is the coordinate of the left face of the piston, and a is the width of the interaction zone. 
Similarly, the gas particle-piston potential for the particles to the right of the piston has the analogous form 

^ ilkfix-Xrf X<Xr, 

[0 X>Xr, ^ ^ 

where Xr = Xrf + a, and Xrf is the position of the piston's right face. We assume that the temperature is low enough 
(or kf is sufficiently large) so that the probability for a particle to reach the piston's surface is negligible. 

In the previous sections it was established that the dynamics of the piston can be deduced from time correlation 
functions describing molecular motion in the field of the piston fixed in space. Below we derive the explicit expressions 
for these correlation functions in the thermodynamic limit, neglecting recollisions of the piston and gas particles due 
to the finite size of the bath. 



6 



Consider the force on the left side of the fixed piston, 

Fi{t) = -kfY^qi{t)0{qi{t)), (32) 

i 

where cji = Xi — Xi is the position of gas particle i relative to the boundary of the interaction zone. 9{x) is the 
step-function, and summation over index i is over all particle in the tube of diameter S to the left of the piston. In 
this section we omit for brevity the subscript for the force on the fixed piston. The simplifying feature of parabolic 
potential is that the time Tc that a particle spends in the interacti on zone of the fixed piston is independent of the 
initial velocity of the particle and given by Tc = tt/w, where uj = sjkjjm. At a given time t, the only gas particles 
in the interaction zone are those that had positive velocities and coordinates q in the interval —vTc < g < at time 
t — Tc- At time the position of the gas particle is determined by qif) = {v/co) smw{t — tin), where tin = t — Tc — q/v 
corresponds to the time at which the gas particle enters the interaction region and q is the position of the gas particle 
at time t — Tc- It then follows that q{t) = (v/lu) sin uj{tc + q/v) = —{v/Lu)smujq/v, which implies that the total 
instantaneous force on the left side of the fixed piston at time t > can be written as 

cx) 

Fi{t) = -kfjdv j dq N{Xi + q,v;t-Tc) ^sin'^. (33) 

-VTc 

In Eq. (33), q = x — Xi and N{x, v; t) is the microscopic linear density of particles defined by 

N{x,v;t) =Y,^{x - Xi{t)') S(v - Vi{t)y (34) 

i 

Similarly, the total instantaneous force acting on the right side of the piston is 

—VTc 

Fr{t) = -kf j dv j dq N{Xr + q,v;t-Tc) ^sin'^, (35) 

-oo 

where q = x — Xj.. 

For a particle outside the interaction zone of the fixed piston (i.e. for x < Xi, and x > Xr), the average linear 
density of particles is {N{x,v)) = nSfM{v), where n is the total (three-dimensional) density of bath particles, and S 
is the cross-sectional area of the piston. It then follows from Eqs. (33) and (35) that the average force acting on the 
left (Fi) and the right (Fr) sides of the piston are (Fi) = —{Fr} = nS/(3. 

It is straightforward to show that the stationary distribution in the vicinity of the fixed piston, including the 
interaction zone, assumes Boltzmann's form 

{N{x, v)) = nSfMiv) eM-U{x)/kBT). (36) 

To calculate the force correlation functions required to analyze the damping terms in the LE, correlation functions of 
the form (A''(Yi)A''(F2) • • • ^{Ys)) must be evaluated, where Y denotes the position-velocity pair {x, v). It is sufficient 
to consider only the case when time arguments are equal for all functions, since time displacement can be transformed 
into spatial displacement for a free particle, i.e. 

N{x,v\t + ti) = N{x -vti,v;t). (37) 

Note that the product NiyifNiy^) can be written as 

N{Y^)N{Y2) =Y.^{yi - Y^)d{Y2 - Y,) 

= J2 ^(^1 - Yi)^iY2 -Yi)+Yl ^(^1 - - Yj). 

i i^j 

Since S{Yi — Yi)S{Y2 — Yi)= d{Yi — 12)^(^1 —Yi), for the ideal gas system one obtains 

{N{Y,)N{Y2)) = 6{Y, - Y2){N{Y^)) + {N{Y,)){N{Y2)). (38) 
For the three-point correlation function, the same arguments lead to the result 
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{N{Y,)N{Y2)N{Ys)) = 5{Yi - Y2)S{Y2 - Y3){N{Yi)) (39) 
+ SiYi - Y2){N{Yi)){N{Ys)) + 6{Yi - Y3) {N (Y^)) {N {Y2)) 
+ S{Y2-Y3){N{Y,)){N{Y2)) + {N{Y,)){NiY2)){NiY3)). 

Eqs. (38) and (39) are the cumulant expansions of {N{Yi)N{Y2)) and {N{Yi)N{Y2)N(Y3)), where the cumulants 

{{NiY,)N{Y2))) = 5{Y, Y2){N{Y,)), 
{{N(Yi)N{Y2)N{Y3))) = S(Yi - Y2)5{Y2 - Y3){N{Yi)), (40) 
{{N{Yi)N{Y2) ■ ■ ■ N{Y,))) = 5{Yi - Y2)5{Y2 -Y^)--- (5(y,_i - Y,){N{Yi)) 

are proportional to the equilibrium density n of gas particles. 

To order A^, the dynamics of the piston is described by the LE (17) with a time-dependent damping coefficient, 

70 (f) = -jI J(| dr(F_Fo(r)), where the evolution of the total force Fo ~ Fi + Fj. is determined by the constrained 
piston-bath Liouville operator Cq. Since {Fi) = —{Fr) and {FiFi{t)) = (FrFr{t)), 

{FFo{t)) = 2{FiFi{t)) - 2{Fif = 2{{FiFi{t))). (41) 

From Eqs. (41) and (33), one sees that (Fi^o(i)) can be expressed in terms of the cumulant {{N{x,v;t — 
Tc)N{x' ,v'; — Tc))). Using property (37) and Eq. (41), the cumulant may be rewritten as 

{{N{x - vt, V] -Tc)N{x', v'; -Tc))) = S{x -vt- x')S{v - v')nSfM{v). (42) 

Then, using Eqs. (41), (33) and (42), one obtains 

{FFo{t)} = —^9{Tc-t)^smujt + Lo{Tc-t)cosujtj J dvfMiv)v^. (43) 
This can be re-expressed in the compact form 

{FFo{t))=Nu,'jio{t), (44) 
where ^o{t) is a dimensionless function (see Fig. 1) given by 

^o(f) = {tc - t) ^sin Lot +{7r- Lot) cos Lot^ (45) 



and 



N = nS^-^ (46) 

CO 



is the average number of particles in the shell of thickness I — a/ (w^) /uj around the piston. The parameter I specifies 
the length at which the average bath particle penetrates the interaction zone, and N = nSl is the average number of 
particles in the layer of thickness I around the piston. Thus, it is evident that in addition to the mass ratio = m/M, 
the relevant physics depends strongly on another characteristic parameter of the system, namely, NX'^, which can be 
interpreted as the ratio of total mass = mnSl of bath particles in the layer of thickness I in the vicinity of the 
piston to the mass M of the piston. 

With these definitions in hand, the time dependent damping coefficient in the LE (17) can be written as 

lo{t) = £ ^ dT{FFo{t)) = uoNCoit), (47) 
and therefore the LE (17) assumes the form 

^^ = XF^{t)-coX^NCo{t)P4t), (48) 
where the damping function (o{t) = oj Jg dr^o(''") is given by 
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(o{t) = \I^^{tc - t)|2(l - cosLot) + (tt - wt) sinwtj + 'i^^it ■ 



) (49) 

(see Fig.l). This expression describes the interesting time-development of the dissipative force, an aspect of the 
dynamics that is outside the scope of more phenomenological models. For t < Tc, Eq. (48) describes essentially 
non-exponential relaxation of the momentum {P{t)) = P(0)e~-^^*^ with 



X{t)=\^Nujj (irCo(T) = A2iVY-|7r + 2wt- (tt- wi)cosu;i-3sina;t|. (50) 

For small Lot, x{t) ^ i^t)"^- On ^ timc-scalc t > Tc the damping function reaches its plateau value and the Markovian 
limit of the Langevin equation, in which the damping coefficient is independent of time, is recovered: 

^Ml = XF\t)-A^^^coX^NP,{t). (51) 

dt V TT 

Note that the characteristic time for relaxation of the momentum Tp = LiJ~^{X'^ N)~^ is governed by the parameter 
A^iV, rather than A^. It is also interesting to observe that the average number of collisions of bath particles with the 
piston for t ^ Tp is of order nS{v'^)^/'^Tp = A"^ and depends neither on temperature, nor on the piston's size, but only 
on the mass ratio. 

Equation (51) is obtained under the condition that the characteristic time Tc for the force autocorrelation function 
is negligible on a time scale Tp of dynamics of the momentum of the piston. Assuming in addition that the random 
force in this equation is Gaussian, one can obtain the Fokker-Planck equation for the momentum distribution function 



dfjP) 

dt 



This coincides exactly with the equation for the piston interacting with the bath particles through a hard-wall potential 
previously obtained by Van Kampen from the master equation [10]. 

The assumption of a Gaussian random force appears to be justified for Eq. (51) describing dynamics on time scales 
much longer than Tc- In this case one can use a coarse-grained description of the dynamics with time resolution 
Tc <C <C Tp. The coarse-graining procedure corresponds to replacing the instantaneous random force in Eq. (51) 
by its average over a time window of duration Af, i.e. F{t) F{t) = At~^ It~^^* ^(*) dt. As previously discussed, 
the number of collisions of bath particles with the piston for the time interval t,, is of order A^^ S> 1. Therefore the 
resolution time interval At may be chosen sufficiently long that the piston experiences many uncorrelated collisions 
during At. Then, according the central limit theorem, one may expect that F{t) is Gaussian-distributed. 

For the more general LE (17) with time-dependent damping, the random force is generally not Gaussian-distributed. 
However, one can easily demonstrate that the distribution of the random force is approximately Gaussian in the limit 
N ^ 1 where the piston interact simultaneously with many bath particles. In fact, the cumulant expansion (23) of 
the multi-time correlation function C2s = {F(ti)F{t2) ■ ■ ■ F(t2s)) contains the products of s pair correlation functions 
{F{ti)F{tj)) ~ {{F{ti)F{tj))). Since a cumulant of any order is proportional N, these terms are of order . The 
other terms in the expansion contain fewer factors of the cumulants and therefore fewer factors of N, and hence may 
be neglected. Then C2s can be approximately expressed as a linear combination of pair correlation functions, a well- 
known characteristic of a Gaussian random variable. The explicit form of the distribution function for the random 
force /(Fq) can be obtained using the inverse Fourier transformation of the generating function (see, for example, 
reference [2]) 

fiPo) = ^ dkexp\-ikFo + f2^-^{{F^))]. (53) 



The cumulants of odd orders vanish, and the cumulant of even orders are ((-Fq*)) = 2{{F^^)) where Fi is given by (33). 
Using (41) one finds that 

where = k'j{v'^) /2(jj'^ , and r(s) is the gamma-function. Substitution (54) in (53) gives the following integral 
representation for the distribution function 
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/(F) = —y rffc exp|-ifcF-\/2^Ar(l-e-(s'=) j| (55) 

If W ^ 1, one can approximately write in the above expression 1 — f^ff'^) « (.9^)^, which leads immediately to 
the Gaussian distribution for the force. The similar arguments hold for the distribution function of higher order 
f{F{h),F{t2),---F{ts)). 



V. THE NONLINEAR LANGEVIN EQUATION 

We now turn our attention to the terms of higher order in A in the A-expansion (13) of the memory function 

K{t) = (FF'^it)) = Y^iKi{t). It was shown in section 3 that for a homogeneous bath Ki{t) = 0, and the first 
non-zero correction to Ko{t) = {FFo{t)) is K2{t) which is of second order in A. From Eqs. (13) and (23) one can see 
that K2 has the structure 

K2{t) = ^A,{t)P; + -A2{t), (56) 
where the functions Ai (t) and A2 (t) are given by 

Ai{t)= [ dh dt2{{GoG2{t,ti,t2))), (57) 
Jo Jo 

A2{t)= f dh r dt2{{GoGo{t-ti)Gi{tM)))- (58) 
Jo Jo 

Substitution oi K Kq + K2 into the exact equation of motion (11) leads to the non-linear generalized LE of the 
form 



dP^{t) 



\F\t)-\^ f dT Mi{T)P.{t - t) - X'^ I dT M2{T)P^{t-T), (59) 

Jo Jo 



dt 

where the memory functions Mi(t) and M2{t) are 

Mi(t) = Mo{t) - -^A,{t) + ^-A2{t), (60) 

M2{t) = -^A,{t). (61) 

It is interesting to note that a nonlinear LE of this form has previously been obtained using a mode coupling ap- 
proach [15]. A similar Markovian version of a nonlinear Langevin equation with a cubic damping term was considered 
by MacDonald [16] on a purely phenomenological grounds. 

Eq. (59) differs from the linear LE (15) with memory fimction Mo(t) = ■^{FFf){t)) not only by the presence 
of nonlinear damping, but also by appearance of correction terms of order A^ to the memory function Mi(t) for 
the linear damping. Note also that the last term on the right-hand side of (59) can be written in the local form 
—\^P^{t) J* dTM2{T) since the non-local correction to this expression A^ J^* dr M2(t) J^_^ dr' P^{t') is of order A^ 
[14]. However, the linear damping term in Eq. (59) cannot be simply written in the local form —\^P^(t) d,TMi{T) 
since in this case the non-local correction has contributions of order A"^ and A* which must be retained. This correction 
can be written in the form 

v2 



/ drMiir) f dr' P,{t') ^ [ drMoir) [ dr'F^'iT') 

Jo Jt-T Jo Jt-T 

-A* f drMoir) f dr' f dr" Mo{r")P,{r' - r") + 0{X^), (62) 
Jo Jt-T Jo 



where we have used the result that Mi{t) = Mo(t) + O(A^) according to (60). The first term in the right-hand side 
of this expression depends on initial coordinates of the bath and may be treated as a small correction to the random 
force F^{t). The second term can be written in the local form 



10 



-X'P^it) [ drMoir) [ dr'f dr" Mo{t") + 0{\^) (63) 
Jo Jt-T Jo 

to order A^. As a result, Eq. (59) can be written in the local form 

^ = AFt(t) - A%i(On(i) - >^'72{t)P!{t), (64) 

with the modified random force 

(i) = +X3 f' drMo (r) /* dr' (r') (65) 

Jo Jt-T 

and the damping functions given by 

li{t)= I dTMi{T)+\^ j drMoir) f dr'f dr"Mo{r"), (66) 
Jo Jo Jt-T Jo 

72(i)= / drM2{r). (67) 
Jo 

For t » rc, the timc-dcpcndcnt coefficients ji{t) attain their limiting timc-indcpcndcnt values 7i, which can be 
obtained from the expressions above by setting the upper integration limit t to infinity. 

It is possible to obtained explicit expressions for the memory functions Mi{t) and the damping functions 7i(t) 
for the extended Rayleigh model. To accomplish this, explicit expressions for the functions Gi defined by Eqs.(19) 
must be computed. It is convenient to express these functions as the sum of two parts corresponding to the force 
acting on left and right sides of the piston, Gi = Gu + Gri- It is somewhat problematic to calculate terms such as 
Gil = e^°'^dFi{t)/dX due to the parametric dependence of Fi{t) on X that is evident when the force is expressed 
in terms of qi{t) and Vi{t) (see Eq. (33)). One straightforward, albeit inelegant, way to circumvent this difficulty is 
to express the force in terms of N{x,v) = N{x,v;t = 0). Details of this technique can be found in the Appendix B. 
Using this approach, we obtain for the left-side part of Gi 



OO 

Gn{ti,t2) = kf j dv j dq N{Xi + q,v;-rc) cos Lo(ti + ^^ 



-v{Tc+tl) 

OO —vti 

— kf9{rc — t2) cosu!t2 J dv J dq N{Xi + q,v;—rc), (68) 

-v{Tc+tl-t2) 

where the integration limit b is 

b = -v{rc + ti- t2)e{rc - t2) - vhe{t2 - n). (69) 
The expression for Gri can be obtained from the one above by replacements, 

/ dv^ dv, / dq^ dq. (70) 

Jo J —OO Ql Q2 

Note that (Gu) = {Gri)- 

For the left-side contribution to G2, we obtain 

00 

Gi2{t,ti,t2) = kf j dv N{Xi-v{t + rc),v;-rc) 


00 6 

■kf Jdv J dq N{Xi + q,v;—rc)— sinuj (t + 

00 

+ kf(j>{tut2) j dv N{Xi-vt,v;-rc), (71) 



-V{t+Tc) 
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where the function (j){t,ti,t2) is 



4>{ti,t2) = 6{t2 - Tc) - 0{ti - Tc)6{Tc - C0SU>t2 

+ 6{tc — ti)6{Tc — t2) cosa;f2 coswti, (72) 



and the upper integration limit b is 

b = -v{Tc + t - t2)0{Tc - t2) - Vte{t2 - Tc). (73) 

A similar expression for Gr2 can be obtained from that for Gn using replacements (70), and, in addition, by 
multiplying the first and the last terms on the right-hand side of Eq. (71) by —1. Recall that we anticipate {G2) = 
by symmetry so that {Gr2) = —{G12), which can be explicitly verified from the expressions above. 

It has been shown above that the damping forces in the nonlinear Langevin equation can be expressed as integrals 
of the cumulants {{G0G2)) and ((GoGoGi)). From Eqs. (68), (71), (33), and (41), one can get an explicit expressions 
for these cumulants. For t > ti > t2 we find 

((GoG2(i,il,i2))) = 2{{GioGi2it,h,t2)}) 



—MtMM) I dvfM{v)v, (74) 
^ Jo 



where 

(j)i{t,ti,t2) = 29{tc — t)smujtcosivticosu)t2- (75) 

The second cumulant required is 

{{GoGoit - h)Gi{t,t2))) = 2{{GioGioit - h)Gn{t,t2))) 

/ k \^ 

= -nSi^-^j Mt,ti,t2) J dvfMivy, 

where the function (j)2 has the form 

^!>2(i, ^1,^2) = 0{Tc — t)cOSLjt2X 

^co{tc — t) cosw{t — ti) + ^ smcij{t + ti) + ^ smu{t — ti)|. 
These results allow to calculate the functions Ai{t) defined by Eqs. (57) and (58), 

Ai{t) = m^io^NUt), (76) 
A2{t) = -mWr'N^2it), (77) 

where dimensionless functions ^,(t) are given by 



a(t) = -^0{To-t)sm^ojt, 

V2lT 



6(i) = -^0{Tc-t)\sm^u)t + u;t{Tr -Lvt)smu}t]. (78) 



/27r 

The memory functions in the non-Markovian LE equation (59), take the form 

M,{t) = Nu;^{(o{t) - 2X^^,{t) - XH2{t)}, (79) 

M2{t)=Nu;^^^i{t), (80) 
m 

Here the function ^o{t) is given by (45) and governs the memory function Mo{t) in the linear LE (15), Mo{t) = 
Noj'^^Qit). The additional A^-corrcction terms leads to a faster decay of the linear damping kernel Mi{t) compared 
to Moit). The kernel for nonlinear damping M2{t) is not decaying function of time but rather has a maximum at 

t = To/2. 



12 



Explicit expressions for damping functions 7i(t) in the local- in-time LE (64) by integrating Mi{t) according to Eqs. 
(66) and (67), and the nonlinear LE for the Rayleigh model with a parabolic potential takes the form 

^ = XF^t) - NX'u:{a{t)P.{t) + \'lc,it)P^it)}, (81) 
where the nonlinear damping function ^2 (i) is given by 



C2{t) =ioJ^ dTii{T) = -^0{Tc - - COSWi + i COS^ Wtj + \\[^0{t - Tc), (82) 

and the linear damping function 71 {t) can be written as 

Ci(t) = Co(i) + A2£i(0 + NX^e2{t). (83) 

The main contribution to the overall damping coefficient Ci(i) is given by the function C,o{t) = co rfr^o(T), while the 
corrections arising at higher orders in the A-expansion are 

£i(t) = -a;^*dr{2a(T)+6(T)}, (84) 

e2{t)=aj' f drUr) f dr' f dr'^r")- (85) 
Jo Jt-T Jo 

Note that ei(t) and £2(4) are of different signs for all t. If N <^ 1 the correction is determined mostly by £i(t) < 
and tends to decrease the damping fiuiction Ci(i)- In contrast, if A'' ^ 1 the main correction comes from e2(i) > 
effectively increasing the linear damping. 

For a coarse-grained description on the time-scale t ^ Tc with a time resolution Tc <C At <C Tp, one can replace the 
damping functions 7,(f) in the LE (81) by their limiting values 7i(rc). 



VI. CONCLUDING REMARKS 



The notion that the character of Brownian motion of a finite-sized particle may depend on parameters other than 
the mass ratio dates back to Lorentz and has been examined by many authors (see [17-19] and references therein). 
It is known that when hydrodynamic effects are important another relevant parameter is the ratio of the mass density 
of the bath to that of the particle. In this paper we have demonstrated that even when hydrodynamic effects are 
absent, as in the extended Rayleigh model, the character of the behavior of a tagged particle may not governed by 
A^ but by the renormalized parameter XI = NX^ , which can be interpreted as the ratio of the average total mass of 
particles in the interaction zone M* to the mass M of the tagged particle. When the average number of particles 
in the interaction zone is large (i.e. A'^ ^ 1), A* ^ 1 is a necessary condition for the applicability of a conventional 
perturbation scheme of derivation of the LE. In this case the conventional assumption of Gaussian random force is 
justified for any time scale. U N < 1, the parameter of the expansion is A^, and the Gaussian force approximation 
holds only on a time scale that is much longer than the characteristic time for the relaxation of the bath. 

Although this paper focuses on the specific model of an ideal gas bath interacting with a Brownian particle through 
a quadratic repulsive potential, many of the results obtained are quite general. In particular, the LE (64) and the 
expressions (66)-(67) for the damping coefficients in terms of microscopic time-correlation functions, which may be 
considered as the generalized version of the fluctuation-dissipation theorem, are limited neither to the specific form 
of interaction potential between the bath and the tagged particle nor to the ideal gas bath. The results of section 
3 concerning the cumulant expansion of the kernel K{t) = {FF^{t)) are also general and not limited to any specific 
model. Combined with quite general theorems about cumulant properties [20], these results may be useful for more 
realistic models with interacting bath particles. 

The explicit expression for the kinetic coefficients and memory functions appearing in the Langevin equations have 
been derived in this paper in the thermodynamic limit, so that any correlations due to finite size of the system are 
neglected. It should be mentioned, however, that the equations themselves, as well as the fluctuation-dissipation 
relations relating the kinetic coefficients to correlation fmictions, also hold for a system with finite baths. For finite 
systems the explicit form of the kinetic coefficients may be rather complicated even for a bath composed of ideal gas 
particles. 

In this paper we have considered the case of the totally symmetric bath when thermodynamic and microscopic 
properties of the gas to the left and to the right of the piston are the same. Some interesting physical implications 
arising in the case of an asymmetric bath will be presented elsewhere [21]. 
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APPENDIX A 

In this Appendix, for the sake of completeness, we show explicitly that the function C4 defined by Eqs. (14) is of 
second order in cumulants of products of G,. C4 can be written as the sum 

where the first constituent, 



{GoGoit - ti)){Go{t - t2)G2{t,t3,U))}, 



is obviously quadratic in cumulants. 



(1) _ 



{{GoGoit -t2))){{Go{t-h)G2{t,h,U))) 
+{{Go{t - h)Go{t - t2))){{GoG2{tMM))) 

+ {{GoGo{t - ti)Go{t - t2)G2{t, t3,t4)))}. 



The second term is 



^[{GoGoit - ti)Ai) - {GoGoit - ti)){Ai) 
-{Giit,t4)) [{GoGoit - ti)Giit -t3,t2- is)) 
-{GoGoit - t,)){G,it - t3, t2 - is))] }, 
where Ai = Sit — ^2)^0(^2 — is)Gi(i2, ^4)- Noting that Ai can be written as 

= Gi(i - is, i2 - i3)Gi(i, i4) + Go(i - i3)G2(i, i2, i4), 

and recalling that due to symmetry (G2(i)) and (GoGi(i)) are zero at all times (see Eq. (21)), one can see that C[ 
is also quadratic in cumulants. 



(2) 



G 



(2) 



^{(Gi(i - i3,i2 - is))((GoGo(i - ii)Gi(i,i4))) 

+ ((GoGo(i-i3)))((Go(i-ii)G2(i,i2,i4))) 
+ ((GoG2(i,i2, i4)))((Go(i - ii)Go(i - is))) 
+((GoGo(i - ii)Gi (i - is, i2 - i3)Gi (i, i4))) 



+ 



The third term is 



where 



((GoGo(i-ii)Go(i-is)G2(i,i2,i4)Gi(i,i4)))}. 
^{(GoBi) - (GoB2)(Gi(i2,i4)) + 2(GoB3)}, 



B, = Sit - t,)Sit ~ i2)Go(i2 - is)Gi(i2, i4), 
B2 = ^(i-ii)^(i-i2)Go(i2-is), 
Bs = Sit - ii)Go(ii - i2)G2(ii,i3, i4), 
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can be expressed as 

Si = G2{t- ts, h - - h)Gi{t, U) + Gi{t- ts, t2 - h)G2{t, h,U) 
+ Gi{t- t3,h - h)G2{t,t2,U) + Gait - h)Gs{tM,t2M), 

B2 = G2{t — ts, tl — ts, t2 — ts), 

B3 = Gl{t - t2,tl - t2)G2it,t3,U) + Go{t - t2)G3it,tl,t3,U). 

(3) 

Then using the symmetry property (21) one can see that C4 is quadratic in cumulants. 



The remaining term 



^|4) ^ J n\ (^GoGi{t,h,t2,t3,ti)) 

' m I 



[{GoGo{t-h)G3it,t2,h,t4)) 
{GoGo{t-h)){G3{t,t2,t3,U))} 



is clearly linear in cumulants, 



= (^^) {{GoG4it,h,t2,h,U))) 

3p2 

+ ^{{GoGo{t - h)Gs{t,t2,t3,t4))). 



APPENDIX B 



In this Appendix the functions Gi defined by Eqs. (19) are evaluated for the; extended Raylcigh model of diffusion. 
These functions arc defined in terms of powers of alternating operators d/dX and mentioned in 

the text, the representation (33), (35) for the force on the fixed piston is unwieldy since it involves qi{t) and Vi{t) 
which arc fmictions of X. For the purpose of evaluating the Gj, it is convenient to express the force in terms of 
N{x, v) = N{x, v;t = 0). For t > 0, we have 



o(t) 



■smoj [t A — 



Flit) = kf Jdv J dq N{Xi + q,v) ^ 

-vt 

00 00 

+e{Tc/2-t)kf j dv j dq NiXt+q,v) q{t) 



00 

+0(re/2 - t)kf J dv J dq N{Xi + q, v) q{t) 



-00 Q(t) 

Q{t) 



+9{t - Tc/2)e{Tc - t)kf j dv j dq N{Xi + q, v) q{t), (Bl) 







where 



a{t) = -ve{t-Tc){t-Tc), 

, , V . 

q{t) = qcoswt H smut, 

CO 

V 

Q{t) = tanut. 

The first term in Eq. (Bl) describes the contribution to the force Fi{t) from the particles that were outside the 
interaction zone at t = (i.e. g < 0) and are in the interaction zone at the moment t (i.e. q{t) > 0). 
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The remaining terms give the contribution from particles that were in the interaction zone at t — (i.e. q > 0) 
and arc still there at the moment t (i.e. q(t) > 0). In fact, all particles in the interaction zone with positive initial 
velocities at t — will be still in the interaction zone a,t t < Tc/2 (the second term), while the particles with negative 
initial velocities will be in the interaction zone at time t < Tc/2 only if at t = they reside deep inside the interaction 
zone, namely q > Q{t) (the third term). For Tc > t > Tc/2 only the particles with positive initial velocity will be in 
the interaction zone at time t provided their initial coordinates are less than Q{t) (the last term). 

Prom expression (Bl), one easily calculates dFi{t)/dX writing dN{Xi + q,v)/dX = dN{Xi+q)/dq and integrating 
by parts to obtain, 

^^ = -kf cos cotNS) (B2) 

ait) 



~kf J dv J dqN{Xi+q,v) — cos lo (t + —"^ . 



-vt 

Here Nz{t) is the number of particles which were in the interaction zone at t = and remain at t > 0, 



N,it) = e{Tc/2 -t)jdvjdq N{Xi + q, v) 





oo 



+e{Tc/2-t) J dv J dqN{Xi + q,v) 

-oo Q(t) 

Q{t) 



+e{t - Tc/2)0{tc -t) j dv J dq N{Xi + q,v). (B3) 





Nz{t) can be written in compact form in terms of the density at time —Tc according to, 

oo —vt 



N^{t) = 9{tc -t) j dv j dq N{Xi + q, v; -Tc). 



(B4) 



—VTc 

In fact, the number of particles in the interaction zone at t = is 

oo 



j dv j dq N{Xi+q,v;-Tc), (B5) 



—VTc 



while at time t it is given by 

oo 

j dv j dq N{Xi + q,v;t-Tc) 



—VTc 



OO —vz 

= J dv J dq N{Xi+q,v;-Tc). (B6) 



+ 

By definition Nz{t) involves the particles which contribute to both integrals (B5) and (B6), which leads to Eq. (B4). 
Combining Eqs. (B2) and (B4), we have 

oo a(t) 



dX 



J dv J dq N{Xi + q, v) cos w (^t + ^) 

-vt 

oo —vt 

kf9{Tc- t)cosujt j dv j dq N{Xi + q,v;-Tc). (B7) 



—VTc 
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For comparison, let us calculate the average derivative of the total force Fo{t) = Fi + Fr, 

dFo{t)\_ /dFi{t)\ 
dX / \ dX / 

_ ^ 2kfnS J ^^^^ dv6{Tc — t) (sin ijt+ {n — tot) cos Lotj . 

Comparing this expression with Eq. (43) for the memory kernel Ko{t) = {FFo{t)), it is evident that {dFo{t)/dX) = 
—/3{ FFo{t)). This is the general result used in the main text, Eq. (10) confirmed using the explicit expression for 

dFo{t)/dX. 

^From Eq. (B7), we obtain the following expression for Gn: 

„£o(ti-t2)^^lM 



dX 

a{t2) 



— J J dq N{Xi + q,v;ti — t2) cosui (t2 + ^ 

-Vt2 

OO —Vt2 

- kfe{Tc~t2)cosu!t2 j dv j dq N{Xi+q,v;ti-t2-Tc), (B8) 

—VTc 

where, according to (13), it is assumed that ti > t2- Expressing this in terms of the microscopic density at time — Tc 
using (37), one obtains Eq. (68) of the main text. 

To evaluate G2, one has take the derivative of Gi with respect to X. Let us express Gi in terms of the density at 
time t = Q, N{x,v), as done above for Fi{t). The first term in the right-hand side of Eq. (B8) involves only particles 
located outside the interaction zone, so using the property (37), the first term can be written as 

00 a{t2) 

— kf Jdv J dq N{Xi+q — v{ti—t2),v) cosw + 

-Vt2 

00 b 

= -kf j dv j dq N{Xi + q,v) cosw {ti + ^ , (B9) 

-vti 

where the upper integration limit is 

h = -V{ti - t2)e{Tc - t2) - V{tl - Tc)e{t2 " Tc) . 

The second term on the right-hand side of Eq. (B8) can be expressed as 

00 —Vt2 

— k}6{Tc — ti) COS ujt2 j dv j dq N{Xi + q — v{ti — t2),v;—Tc) 

—VTc 

00 —Vtl 

= -kfe{Tc - ti) coswti j dv j dq N{Xi + q,v; -Tc). (BIO) 

-v{Tc+tl-t2) 

If ti > Tc- then all particles which contribute this integral at t = will be outside the interaction zone in the g-interval 
from ~v{ti — ti) to —v{ti — Tc). Therefore, for ti > Tc the second term equals 

00 -V(ti-Ta) 

-kf0{Tc-t2)9{ti-Tc)cosLot2 J dv J dqN{Xi + q,v). (Bll) 

-v(ti-t2) 
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If ti < Tc, then two sets of particles contribute to expression (BIO). The first set of particles is composed of particles 
which at time —Tc are in the g— interval from —v{tc + ti— ^2) to —vTc- At t = 0, these particles will be outside the 
interaction zone in the interval {—v{ti — 12), 0), and hence their contribution is 

00 

-kf cos LJt2 j dv j dqN{Xi+q,v). (B12) 

-v{tl-t2) 

The second group of particles arc those that were in the interval from —vt^ to —vti at time ^r^. At i = 0, all 
these particles will be in the interaction zone. Taking into account Eq. (B4), the corresponding contribution can be 
written as —kf cosu)t2Nz{ti). Using expression (B3) for Nz{t), one arrives at the following representation for Gn in 
terms of N{x, v), 

oc b 

Gn{ti,t2) = -kf j dv j dq N{Xi+q,v) cosLj(ti + ^^ 

-vti 

00 -f(tl-Tc) 

—kf6{Tc — t2)0{ti — Tc) cosu!t2 Jdv J dqN{Xi + q,v) 

-v{ti-t2) 

00 

—kf6{Tc — t2)0{Tc — h) cosu!t2 J dv J dqN{Xi+q,v) 

-v{ti-t2) 

00 00 

-kfe{Tc-t2)0{Tc/2-t-i)cosut2 J J dq N{Xi + q,v) 



00 

-k f e{Tc- 12)6 {Tc/2-ti) cos wt2 j dv j dq N{Xi+q,v) 



-00 Q(ti) 

-kf6{Tc - t2)0{Tc - tl)6{tl - Tc/2) C0SU)t2 X 

00 Q(ti 





Taking derivative of this expression with respect to X gives 



j dv j dqN{Xi + q,v). (B13) 



00 

g^Gnih,t2) = kf I dv N{Xi - vh,v) (B14) 


00 

+kf[9it2 - Tc) - e{Tc - t2)e{ti - Tc) coswta] j dv N{Xi - v{ti - Tc), v) 



00 b 

-kf J dv J dq N{Xi +q,v) ^ sinw (ti + ^) 

-vti 



+e{Tc/2 - ti)e{Tc - t2)kcosujt2 j dv N{Xi + Q{ti),v) 

—00 

00 

-e{Tc-ti)6{ti-Tc/2)e{Tc-t2)kcoswt2 j dv N{Xi + Q{ti),v). 
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The last two terms can be written as 



CXD 





oo 

Vt-] V 



sinwfi ' sinwfi ' 



oo 

= kf6{Tc — t2)6{Tc — ti) coswticosut2 J dv N {Xi — vti, V, —Tc) , (B15) 



where we have used the property that if the initial coordinate and velocity of the particle in the interaction zone are 

q{0) = -, v{0) = --^ (B16) 
u> tanojt 

with V > and < t < Tc, then a,t t = —Tc 

q{-Tc) = -v{-Tc)t, v{-Tc) = (B17) 

SVSlLOt 

Substituting (B15) into (B15), acting on the result by the propagator e^f(t-*i)^ an^i using again the property 
N{x, v;t + T) = N{x — VT, v; t) for the motion outside the interaction zone, we finally obtain for t > ti > t2, 

Gz2(t,ii,i2) = e^°(*-*^)^Gi(ti,t2) 



= kf J dv N{Xi - vt, v) (B18) 



oo 

+kf0{h - Tc)9{t2 - Tc) j dv N{Xi - V{t - Tc), V) 


oo 

— kf9{Tc — t2)0{t\ — Tc) C0SL0t2 j dv N{Xi — v{t — Tc), v) 



oo h 

-kf j dv j dq N{Xi +q,v)'^ sin LO (t+^^ 

-vt 

oo 

+kfO{Tc — t2)0{Tc — ti)coswti cosa;f2 J dv N{Xi — vt,v;—Tc), 



where b = —v{t — t2)0{Tc — — v{t — Tc)9{t2 — Tc). Expressing the first three terms through A''(a;, v; —Tc) rather than 
N(x,v), one arrives at Eq. (71) of the main text. 
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FIG. 1. The function £,o{t) (solid line) governing the time dependence of the memory function Mo{t) = Ncj^S^oit) in 
non-Markovian Laiigcvin equation (15), and the damping function Co(^) (dashed line) in the local Langevin equation (48). 
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